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The best way of controlling the theoretical uncertainties of the perturbative QCD predic- 
tions is the direct analytical or numerical calculation of the concrete terms in the corresponding 
perturbative series. However, after the results of the calculations of the next-to-next-to-leading 
order (NNLO) corrections to the number of physical quantities became available fip^d]) ex- 
perimentalists and theoreticians started to be interested in the effects of still uncalculated 
higher-order terms. In the work of Ref. f^, two "optimal" methods of fixing the renormal- 
ization scheme ambiguities were used to estimate the next-after-next-to-next-to-leading order 
(NANNLO) corrections to certain renormalizat ion-group invariant quantities. These methods 
were the principle of minimal sensitivity (PMS) and the effective charges approach (ECH) 
0, which is known to be a posteriori identical to the so-called scheme- invariant perturbation 
theory The quantities studied in Ref. are the e"'"e~-annihilation ratio R{s), the r-lepton 
decay ratio Rr, and the Bjorken non-polarized and polarized sum rules. 

However, the quantities obeying the renormalizat ion group equations with anomalous di- 
mension functions were not considered in Refs. [0, 0. In this note we will fill in this gap and 
apply the generalization of the ideas used in Ref. ||^ to estimate the higher-order corrections 
to the singlet part of the Ellis- Jaffe sum rule (EJSR), recently calculated at the O(a^) order 
0. We present the concrete O(a^) estimates for / = 3 numbers of flavours in two related 
forms, namely in the factorization-scheme-invariant form and in the form that necessitates the 
application of an additional guess about the value of the unknown four-loop coefficient of the 
corresponding singlet anomalous dimension function. 

The experimental measurements |10] of the structure functions g^^"^ of the polarized deep- 
inelastic lepton-nucleon scattering has stimulated number of works aimed at a theoretical study 
of the structure function g^^^^ and their first moment, namely the EJSR (see e.g., ||11||-||15||). 



The theoretical expression for the EJSR can be presented in the following form: 

EJSRiQ') = l\l^''\x,Q^)dx = EJ^sm + EJsim . (1) 

JO 

The non-singlet contribution to this sum rule is defined as 

EJNsm = (1 - « - E d^^'^') (±^«3 + 4^8) , (2) 

^ 12 36 

where a = Os/tt, = Au — Ad, as = Au + Ad — 2 As and Au, Ad, As can be interpreted as 
the measure of the polarization of quarks in a nucleon. 



The singlet contribution to Eq. (|T|) reads 



/a(Q^) 'yr.r(r) 1 
^^rfx}-ASi„„ , (3) 



where 



I3{x 



AJ^inv = exp{- J ^^^^^(ia;}AS(/i^^ 



1 



and AS = A-u + Ad + As. The first three coefficients of the QCD /5-function 



/3(a) = -5:Aa^^^ (4) 

i>0 

and of the anomalous dimension of the singlet axial current 

lsi{a) =Y.l^d^' (5) 

j>0 

are known in the MS scheme from the results of calculations [T^, [|17| respectively. They have 
the following numerical form: 

(3o = 2.75 - 0.167 / 

= 6.375 - 0.792 / 
p2 = 22.320 - 4.369 / + 0.094 f (6) 

and 

70 = 

71 = -0.5/ 

72 = -2.458 / + 0.028 f . (7) 

The coefficient 73 is unknown and we will need to somehow fix its value in the process of further 
considerations. 

The perturbative expression for the coefficient function 

C{a) = l + J2r,a' (8) 



is explicitly known in the MS scheme at the O(a^) level from the results of Refs. [|l8l. The 
results of direct calculations are 

ri = —1 

r2 = -4.583 + 1.162 /. (9) 

Our aim is to estimate the value of the coefficient r^, in the MS scheme , thus roughly fixing the 
uncertainties due to the lack of knowledge of the explicitly unknown correction to the singlet 
part of the EJSR at the 0(0;^) level, already achieved in the direct calculations of the related 
non-singlet contribution This estimate will be obtained with the help of a procedure that is 
analogous to the one of Ref. which was used in Ref. to estimate the O(a^) contribution 
to the non-singlet part of the EJSR. As was already demonstrated in Ref. this procedure is 
working quite well at the 0{al)-\eYe\. Indeed, the obtained in Ref. estimates of the order 
0(Q;f)-coefficients of the quantities under consideration are in good agreement with the results 
of direct calculations of Refs. This fact stimulates further applications of the methods 

used. 
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Let us follow the considerations of Ref. [T^ and define the renormalization-group-invariant 
quantity: 



where 



and 



R 



EJsi 



"0 



Q 



2 d\n[EJsim] 



75/ (a) + P{a) 



dC{a)/da 



C{a) 



-2/5or2 + /3or? - /^in + 72 
-3/5or3 + 73 - /(/3o, /32, n, r2' 



(10) 



(11) 



/(/3o, A, /32, ri, r2) = /3or? - 3/3orir2 + 2Ar2 - ^r^ + (3^r^ . (12) 

Using the ideas of Refs. 0, pO[, 0, one can insert the values of the coefficients dQ,di and 
ci = /9i//?o into the following expression: 



SIjSI 



d^'d 



2 -d','df\-df + c,). 



(13) 



This can be obtained as the residual term in the re-expansion of the effective charge Rejsi = 
^o^^'ech in terms of the coupling constant of the initial scheme, in our case the MS scheme. The 
result of re-expansion of a similar expression obtained within the framework of the PMS ap- 
proach will differ from Eq. (jl^) only slightly, namely {dQ^d^^)pMs = {dQ^d2^)ECH + {'^/^)dQ^c\. 
Therefore, we consider the term of Eq. (13) as the one that simulates the main contribution of 
the non-calculated NNLO correction. 



As the result of apphcation of this procedure we get the following expression of the coefficient 



d^ (ig 



+ 



73 /(/3o,/5i,/32,ri,r2) 



(14) 



3/3o 3/5o Wo 

where the anomalous-dimension term 73/(3/3o) remains unknown. Note that this unknown term 
cancels in the expression for Eq. (^ after taking into account the corresponding expansion of 
the anomalous- dimension term 



exp 



dx 



1 



7i 



+ 



+ 



7? _ 72 7iA 
PI A 




73 _ 71/^2 ^ 1if3f _ 72/^1 
Pi Pi Po , 



a 



(15) 
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Taking now / = 3 numbers of flavours and using Eqs. (j^)- (|) and Eqs. (|TT])-(|T^) we obtain 
the following expression for the singlet contribution to the EJSR: 



1 - 0.333a - 0.549a^ - 2a^ 
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(16) 



where the 0{a) and O(a^) corrections are known from the results of explicit calculations (see 
E.g., Ref. IP) and the O(a^) contribution is our estimate of the value of a correction that is still 
not calculated. Notice once more, that the possibility to apply the above-described approach for 
the estimates of the O(a^) corrections to the singlet part of the EJSR is supported by the good 
agreement of the results of similar estimates of the 0{al) corrections to the corresponding non- 
singlet contribution for / = 3 numbers of flavours with the results of the explicit calculations 
of Ref. §. 

In order to get a similar estimate for the Q^-dependent normalization of AE, namely for the 
case when AS = AS(/x^ = Q^), it is necessary to somehow flx the value of the unknown 73 term 
of the singlet anomalous dimension. For / = 3, we will use the following bold guess-estimate: 



73 



71 



-34, 



(17) 



which is supported by the detailed consideration [gl[] of the results of the NNLO calculation of 
the non-singlet anomalous dimension for the flrst four even moments ||22|| . 



Using Eq. ([TtD we obtain the bold guess- motivated estimate of the 0{a^) contribution to 
the singlet part of the EJSR in the case when the dependence of AS is specifled: 



EJsiiQ' 



1.096a^-3.7a'' 
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AS(g2 



Note, however, that we are unable to present a similar estimate of the O(a^) correction to the 
singlet part of the EJSR for / = 4 numbers of flavours. Indeed, for the case of / = 4 the 
value of the coefficient in Eqs. ([10|) and ([Tl| ) is almost nullified (rfg^ ~ since ri ^ 71 //5o). 
Therefore, it is impossible to determine the value of the corresonding coefficient df^ from the 
expression for d^^df^ presented in Eq. (|ll]). This example demonstrates the limitations of the 
procedure discussed above. 

To conclude, we touched the problem of fixing the uncertainties due to still-uncalculated 
0{al) corrections to the singlet contribution into the EJSR, for / = 3. Combining the obtained 
estimates with the results of available NNLO calculation Q of the non-singlet contributions 
into the EJSR and with the corresponding NANNLO estimates |^, we arrive at the following 
expressions for the EJSR in the MS scheme , related to / = 3: 



g{^''\x,Q'')dx 



+ 



1 - a - 3.583a2 - 20.2150^ - 130a^ 

1 



1 - 0.333a - 0.549a^ - 2a'' 
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AS- 



(19) 
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or 



Jo 



1 - a - 3.583a2 - 20.215a^ - 130a^ 



, 1 1 
* 12 36 



+ 



l-a- 1.096a^-3.7a'' 



9 mQ') 



(20) 



It can be seen that the perturbative contributions to the singlet part of the EJSR, includ- 
ing the 0{al) term that we estimated, are negative. They are significantly smaller than the 
coefficients of the perturbative series of the non-singlet part, which include the results of the 
concrete 0{al) calculations |Q and the estimates of the O(a^) terms 0, that are in qualitative 
agreement with the results of application of the Pade resummation technique p3|. The concrete 
physical applications of the results obtained are discussed in Ref. ITH . 
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